Quantum q-breathers in a finite Bose-Hubbard chain: The case of two interacting 

bosons 
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We study the spectrum and eigenstates of the quantum discrete Bose-Hubbard Hamiltonian in a 
finite one- dimensional lattice containing two bosons. The interaction between the bosons leads to an 
algebraic localization of the modified extended states in the normal mode space of the noninteracting 
system. Weight functions of the eigenstates in the space of normal modes are computed by using 
numerical diagonalization and perturbation theory. We find that staggered states do not compactify 
in the dilute limit for large chains. 
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I. INTRODUCTION 

The study of discrete breathers in different physical 
systems has had remarkable developments during the last 
two decades [2, H, H, El- These excitations are generic 
time-periodic and spatially localized solutions of the un- 
derlying classical Hamiltonian lattice with translational 
invariance. Their spatial profiles localize exponentially 
for short-range interaction. Recent experimental obser- 
vations of breathers in various systems include such dif- 
ferent cases as bond excitations in molecules, lattice vi- 
brations and spin excitations in solids, electronic cur- 
rents in coupled Josephson junctions, light propagation 
in interacting optical waveguides, cantilever vibrations 
in micromechanical arrays, cold atom dynamics in Bose- 
Einstein condensates loaded on optical lattices, among 
others 

H 1 S 1 1 M El HI El • In many cases quan- 
tum dynamics is important. Quantum breathers con- 
sist of superpositions of nearly degenerate many-quanta 
bound states, with very l ong times to tunnel from one 
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lattice site to another 

m, m h n m, m" 



breathers, though being extended states in a translation- 
ally invariant system, are characterized by exponentially 
localized weight functions, in full analogy to their classi- 
cal counterparts. 

Recently the application of these ideas to normal 
mode space allowed to explain many facettes of the 
Fermi-Pasta-Ulam (FPU) paradox (3(| , which consists of 
the nonequipartition of energy among the linear normal 
modes in a nonlinear chain. There the energy is localized 
around the initial normal mode which is excited. Intro- 
ducing the notion of q-breathers [3J H3, HH , which are 
time-periodic excitations localized in the normal mode 
space, the FPU paradox and some related problems were 
successfully explained. Despite the fact that the interac- 
tion in normal mode space is long-ranged, it is selective 
and purely nonlinear, thus q-breathers localize exponen- 
tially in normal mode space. 



In this paper we address the properties of quantum q- 
breathers. We study a one-dimensional quantum lattice 
problem with two quanta. By defining an appropriate 
weight function in the normal mode space we explore 
the localization properties of the eigenstates of the sys- 
tem. We observe localization of the weight function as 
a function of the wave number, which we interprete as 
a signature of quantum q-breather excitations. By using 
a numerical diagonalization of the Hamiltonian and non- 
degenerate perturbation theory we find algebraic decay 
of the weight function in the normal mode space, at vari- 
ance to the exponential decay found for q-breathers in 
the case of a classical nonlinear system. Another intrigu- 
ing difference is based on the interference effects of two 
interacting quanta. For the general case the quantum q- 
breather states approach the noninteracting eigenstates 
in the dilute limit of large chains. However, states with 
Bloch momentum close to ±7r keep their finite localiza- 
tion in that limit. 

In section II we describe the model and introduce the 
basis we use to write down the Hamiltonian matrix. In 
section III we review results on the properties of two- 
quanta bound states - the simplest versions of a quantum 
breather. In section IV we consider the case of extended 
states. We introduce a weight function to describe local- 
ization in the normal mode space, and obtain analytical 
results using perturbation theory. We present our numer- 
ical results obtained by diagonalization of the Hamilto- 
nian matrix, comparing them to analytical estimations. 
We conclude in section V. 



II. THE MODEL 

We study a one-dimensional periodic lattice with / 
sites described by the Bose-Hubbard (BH) model. This is 
a quantum version of the discrete nonlinear Schrodinger 
equation, which has been used to describe a great variety 
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of systems [34| . The BH Hamiltonian is given by [3£ 



structure as in [35| for the case of the BH system: 



where 



H = H + 7.H1 



Ho = -'£ l bf(b j -i + b j+1 ), 
3=1 
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with g = 1 + t and p = T~(f +1 ^ 2 +r~(f~ 1 ^ 2 . By varying 
the Bloch wavenumber in its irreducible range, we obtain 
the eigenenergy spectrum shown in Fig[T] 



Here bj and bj are the bosonic creation and annihila- 
tion operators which satisfy the commutation relations 
[bi,bj] = 8ij, [bf,bf] = [bi,bj] =0. 7 is the parameter 
controlling the strength of the interaction, and the chain 
of length / is subject to periodic boundary conditions. 
The chain is translational invariant and the Hamiltonian 
(HJ) commutes with the number operator N — Y^j=i bfbj, 
whose eigenvalue is denoted by n. We consider the sim- 
plest non-trivial case of n — 2. It is of direct relevance 
to studies and observations of bound two-vibron states 

[H M M HI EE EH, S3, EE El EE EE S3 

In order to describe the quantum states, we use a num- 
ber state basis [IE] \&n) — \ni,n 2 , ■■■,«/), where rii rep- 
resents the number of bosons at site i (n = Y)rii). As 
an example 1 0200000) corresponds to a state with two 
bosons on the second site and zero bosons elsewhere. For 
a given number of bosons each eigenstate is a linear com- 
bination of number states with fixed n. In addition to 
the number of quanta n there are n — 1 further quantum 
numbers which define the relative distance between the 
bosons. For n = 2 that reduces to defining one further 
relative distance j — 1 between the two quanta, which can 
take (/ + l)/2 different values in our case: 



1*2 



3=1 



(4) 



Due to translation invariance the eigenstates of H are 
also eigenstates of the translational operator T, where 
t = cxp(i/c) is its eigenvalue with k = 2m// f being 
the Bloch wave number and v e [— (/ — l)/2, (/ — l)/2]. 
Due to periodic boundary conditions T\m,n-2, ■ ■ ■ , n/) = 
|n/, 774, 77,2, • • • , n/_x). For the sake of simplicity we deal 
with an odd number of sites /. Thus we can construct 
number states which are also Bloch states: 



l^2> = 



|10_^_01 

3-1 



(5) 



With this basis we can derive the eigenenergies for each 
given Bloch wave number k from Hk\^n) = E\^f n ) after 
computing the eigenvalues of the matrix with the same 



III. BOUND STATES: LOCALIZATION IN 
REAL SPACE 

In Fig[TJJa-c) we show that as the interaction parameter 
is increasing, an isolated ground state eigenvalue E 2 (k) 
appears for each k that corresponds to a bound state 
[35| . For this isolated ground state there is a high proba- 
bility of finding two quanta on the same site. In the limit 
/ — > 00 the bound state eigenvalue has the analytical ex- 
pression [H, El| : 



E 2 (k) = -V7 2 + 16cos 2 fc/2, 



(7) 



and the corresponding (unnormalized) eigenvector v = 
(vi,V2, . ..) is (Hi 



v = I -^=,/i,/i 2 ,/i 3 , . . . j . 



where 



(7 + E 2 {k))e lk l 2 
4cos(fc/2) 



(8) 



(9) 



A suitable weight function of this isolated ground state 
has the form: 



C, 
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2(3-1) _ p 2A(j-l) 



3 > 1. 



(10) 



where C\ = 1/2 and A = In Since |/i| 2 < 1 for 7 =^ 0, 
the weight function shows exponential decay when the 
distance between the two bosons increases. That result 
corresponds to the exponential localization of classical 
discrete breathers [E 0, IE 0|- However note that for 
|fc| — ► 7r we have fj, — ► independently on the value of 
7^0. Thus one obtains compact localization. Note 
that it is said that a state is compact in a certain basis, 
if it occupies a certain subspace, but has exactly zero 
overlap with the rest. 

The compact localization for \k\ — * ir is not observed 
in the classical limit, and relies on the fact that the 
Schrodinger equation is a linear wave equation which ad- 
mits (destructive) interference effects. 



3 



(a) 




-1.0 -0.5 0.0 0.5 1.0 
k/7t 



(b) 




-1.0 -0.5 0.0 0.5 1.0 
k/jc 



(C) 




-10 



-1.0 -0.5 0.0 0.5 1.0 
k/it 

FIG. 1: Energy spectrum of the Bose-Hubbard model for dif- 
ferent values of the interaction 7: (a) 7 = 0.1, (b) 7 = 1.0, 
and (c) 7 = 10. Here / = 101. 



IV. QUANTUM Q-BREATHERS: 
LOCALIZATION IN NORMAL MODE SPACE 

All the other states (except the bound state) form the 
two quanta continuum. Their energies for 7 = cor- 
respond to the sum of two single particle energies with 
the constraint that the sum of their momenta equals the 
Bloch momentum k. One arrives at 

E° kM = -2[cos(fc 1 ) + cos(*i + k)], (11) 

where ki — nvi/Kf + l)/2] — k/2 is the conjugated 
momentum of the relative coordinate (distance) of both 
quanta and v\ = 1, . . . , (/+l)/2. E k k has a finite spread 
at fixed k (see FigJT]) . However for k = ±ir the spectrum 



becomes degenerate. Thus for |fe ± tt| <C 1 the eigenen- 
ergies are very close (almost degenerate). Remarkably 
the bounds of the spectrum for 7 7^ are very well de- 
scribed by the 7 = result. Increasing 7 at fixed k, the 
eigenenergies will slightly move, but never cross. Thus 
a continuation of an eigenstate at 7 = to 7 7^ will 
preserve its relative ordering with respect to the other 
eigenenergies. 

For 7 7^ these quantum q-breather states will be de- 
formed. In analogy to the study of the fate of normal 
modes in classical nonlinear systems [3ll l32t |33j , we will 
study the changes of the two-quanta continuum. For fi- 
nite / and 7 the new states will be spread in the basis of 
the 7 = continuum. For / — * 00 one expects that the 
new states become again identical with the 7 = states, 
since the two quanta will meet on the lattice with less 
probability as / increases. Thus we will test the com- 
pactification of the new states in the 7 = eigenstate 
basis both for 7 — * and for / — * 00. 

We compute the weight functions in normal mode 
space in order to probe the signature of quantum q- 
breathers. For this purpose we start by using pertur- 
bation theory to set up these weight functions, where 
Hi is the perturbation. We fix the Bloch momentum k, 
and choose an eigenstate |\E'~ ) of the unperturbed case 
7 = 0. Upon increase of 7 it becomes a new eigenstate 
l^fc ), which will have overlap with several eigenstates 
of the 7 = case. We expand the eigenfunction of the 
perturbed system to the first order approximation: 

The perturbation of strength 7 is local in the matrix rep- 
resentation @ , thus the relevant perturbation parameter 
is 7//. This has to be compared to the typical spacing 
of unperturbed eigenenergies. For Bloch wave numbers 
far from ±7r the spacing is of order 1//, so the approxi- 
mation should work for 7 < 1. For Bloch wave numbers 
close to ±7r the approximation breaks down if 7 > n— \k\. 

The off-diagonal (fci 7^ ki) weight function at the first 
order is given by : 

|(*S IJTiltf? )| 2 

c(k i;kl ) = K^J^)! 2 = ■ (i3) 

1 % 1 fei 1 

and E% are the eigenenergies of the unperturbed 

system. With A = ki — k% the weight function can be 
rewritten in the following form 
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2„,2 



A 2 7 



64(/ + l) 2 cos 2 (|) sin 2 (§) sin(2^) cos(|) + cos(2^) sin(f ) 
where A is a constant. For 7 = 0, l*^) = I*- ), and the weight function is compact. For |A| <C 1 



(14) 



C(fci;fci 



(/ + l) 2 f cos 2 (|)[sin(2ktfc; 



A cos (2fc^) 



(15) 



From this formula we obtain several interesting results. 
First of all, the decay of the weight function with increas- 
ing A means that we have localization in normal mode 
space. For 2k\ + k ^ 0, 2ir, we have 



C 



7 



Cf + 1) 2 A 2 ' 



(16) 



We find algebraic decay ~ 1/A 2 of the weight function, 
and for 7 — » or / — > 00 the weight function compacti- 
fies. For 2ki + k = 0, 27r, we have 



C 



7 



1 



(/ + l) 2 A 4 



(17) 



Here we find algebraic decay ~ 1/ A 4 of the weight func- 
tion, that also compactifies when 7 — > or / — > 00. 
Finally, for k close to ±tt and large /, C ~ 7 2 /A 2 . Thus 
we find that the /-dependence drops out for staggered 
states \k ± vr| <C 1, and these states do not compactify 
for / — > 00. That is a remarkable quantum interference 
property, since both simple intuition (see above) and clas- 
sical theory predict the opposite. 

In Fig[2] we show numerical results obtained by diag- 
onalization of the Hamiltonian for different values of 7. 
We find localization in normal mode space, which can be 
interpreted as a quantum q-breather. When increasing 7 
the quantum q-breather becomes less localized, and for 
large values of the interaction (from 7 = 10 on) results 
do not change. The dashed lines are the results using the 
formula (JT5J) with A 2 = 3.8, value that was obtained by 
fitting the numerical results for the lowest 7 (= 0.001). 
We can see good agreement with numerical results up to 
7 = 1, beyond which perturbation theory does not fit 
anymore. In Fig[3l we show that the weight function is 
more localized for k = and less localized for k — ► — it. 

While probing the influence of the size of the nonlinear 
quantum lattice on the localization phenomenon, we find 
in Fig|4]that as the size increases the states compactify 
as we expected. In Figj5]we see the 1/A 2 decay for eigen- 
states fulfilling 2ki + k ^ 0, 2tt (k = 0), and in FigUthe 
1/A 4 decay for eigenstates fulfilling 2ki + k — 0, 27r (also 
k = 0). Both results agree with the analytical results 
using perturbation theory. 

In FigJT] we observe the predicted independence of the 
localization phenomenon from the size of the system 
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FIG. 2: Weight function for different values of the interaction 
using the formula (|14[1 with A 2 



7. Here f=101, k — 0, and ki — |tt. Dashed lines are results 
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lines are results using the formula (|14[1 with A 2 



7 7r. Dashed 



when k is close to — tt. It is interesting that in this case 
the weight function does not compactify in the dilute 
limit / — > 00 as one would expect from simple grounds. 
The reason is that the larger /, the closer we can tune 
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V. CONCLUSIONS 

In this work we studied the properties of quantum 
q-breathers in a one-dimensional chain containing two 
quanta modeled by the Bose-Hubbard Hamiltonian. To 
explore localization phenomena in this system we com- 
puted appropriate weight functions of the eigenstates in 
the normal mode space using both perturbation theory 
and numerical diagonalization. We observe localization 
of these weight functions, that is interpreted as a signa- 
ture of quantum q-breathers. The localization is stronger 
when the size of the system increases. Unlike the clas- 
sical case where the localization is exponential, here we 
found algebraic localization. This is a long range behav- 
ior, which follows from the fact that the interaction 7 
induces a linear perturbation of the eigenstates which is 
local in real space, and also local in the matrix represen- 
tation in ((5]) • That induces a mean-field type interaction 
between the normal modes, and naturally leads to alge- 
braic localization. Note that the matrix ((SJ) is formally 
analogous to a semi-infinite tight-binding chain with a 
defect at one end. Nevertheless it appears in our con- 
text when starting with a translationally invariant sys- 
tem, but with many-particle states which include inter- 
action. 

Since the effective interaction strength drops in the di- 
lute limit of large chains, we observe stronger localization 
(except for the case of staggered states). The crucial 
difference to the classical model is, that while the lin- 
ear classical dynamics coincides with the single particle 
quantum problem, nonlincarity in the classical model ef- 
fectively deforms the single particle dynamics (and adds 
many other features like chaos etc). The interaction in 
the quantum problem takes the wave function into the 
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new Hilbert space of many-body wave functions, which 
is still a linear space, but higher dimensional. Another 
feature of that quantum interaction is the fact that stag- 
gered states do not compactify in the dilute limit of large 
chains. That property is based on the interference of 
quantum states, and is not observed in the correspond- 
ing classical nonlinear equation. A similar (yet weaker) 
signature of quantum interference is the observed change 
of the power of the algebraic decay from two (generic) to 
four when choosing particular values of the wave number 
k±, which depend on the Bloch wave number k. And yet 
another signature of quantum q-breathers is the fact that 
they keep a finite localization in the limit 7 — * 00 as seen 
in FiglH and at variance to their classical counterparts, 
which turn from exponentially localized to completely de- 
localized in that limit. The reason is that in this strong 
interaction limit extended two-boson states correspond 
to their noninteracting counterparts which are projected 
onto the basis space which does not contain doubly occu- 
pied chain sites, while strong nonlinearity in the classical 
problem completely deforms periodic orbits of the non- 
interacting system. 

We are aware of the fact that the quantum problem 
studied here is a linear one (in terms of differential equa- 
tions). Its correspondence to a classical nonlinear system 



can be observed in the limit of many bosons when treat- 
ing the many particle quantum states within a Hartree 
approximation, which projects onto product states. Of- 
ten the classical description is also achieved using suitable 
(e.g. coherent state) representations. The presented re- 
sults have an unambigous meaning in the chosen basis of 
the noninteracting system. Yet they will of course in gen- 
eral depend on the chosen basis. Therefore it remains a 
puzzling question, how to restore exponential localization 
of classical q-breathers from the algebraic decay of quan- 
tum q-breathers with two bosons, in the limit of larger 
numbers of bosons. The fate of quantum q-breathers in 
higher dimensional lattices is another interesting open 
question, which will be left to future work. 
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